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Abstract 

Multivariate random fields whose distributions are invariant under operator- 
scalings in both time-domain and state space are studied. Such random fields 
are called operator-self-similar random fields and their scaling operators are 
characterized. Two classes of operator-self-similar stable random fields X = 
{X(t),t 6 M. d } with values in R m are constructed by utilizing homogeneous 
functions and stochastic integral representations. 

Keywords: Random fields, operator-self-similarity, anisotropy, Gaussian ran- 
dom fields, stable random fields, stochastic integral representation. 

1. Introduction 

A self-similar process X = {X(t),t € M} is a stochastic process whose finite- 
dimensional distributions are invariant under suitable scaling of the time-variable 
t and the corresponding X(t) in the state space. It was first studied rigorously 
by Lamperti [16] under the name "semi-stable" process. Recall that an M m -valued 
process X is called self-similar if it is stochastically continuous (i.e. continuous in 
probability at each tGR) and for every constant r > 0, there exist a positive number 
b(r) and a vector a(r) € K m such that 

{X(rt),t £l} = {b(r)X(t) + a(r),t G K}, (1.1) 

where = means equality of all finite-dimensional distributions. Lamperti [16] showed 
that if X is proper (see below for the definition) then b(r) = r H for some H > 0, 
which is called the self-similarity index or the Hurst index in the literature. 

Self-similar processes have been under extensive investigations during the past 
four decades due to their theoretical importance (e.g. they often arise in functional 
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limit theorems) and their applications as stochastic models in a wide range of scien- 
tific areas including physics, engineering, biology, insurance risk theory, economics, 
mathematical finance, just to mention a few. 

The notion of self-similarity has been extended in two ways. The first extension 
is to allow scaling in the state space M m by linear operators (namely, b(r) in (II. ip is 
allowed to be a linear operator on M m ) and the corresponding processes are called 
operator-self-similar processes in the literature. More specifically, Laha and Rohatgi 
|15| first extended Lamperti's notion of self-similarity by allowing b{r) in (jl.ip to be 
in the set of nonsingular positive-definite self-adjoint linear operators on M. m . Hudson 
and Mason [12] subsequently allowed b(r) to be an arbitrary linear operator on M. m . 
The operator-self-similarity defined by Sato [28] has an additional assumption that 
a(r) = in (11.11) . Thus the operator-self-similarity in the sense of Sato [28] is 
stronger than that in Hudson and Mason [12]. Various examples of operator-self- 
similar Gaussian and non-Gaussian processes have been constructed and studied by 
Hudson and Mason [12], Sato [28], Maejima and Mason [17], Mason and Xiao [19] . 
Didier and Pipiras [9]. The aforementioned extensions to operator-self-similarity 
is useful for establishing functional limit theorems for multivariate time series and 
their statistical inference [21] . 

The second extension is for random fields (i.e., multi-parameter stochastic pro- 
cesses) which is to allow scaling by linear operators on the multiparameter "time" - 
variable t G M. d . This was done by Bierme, Meerschaert and Schefner [3]. In their ter- 
minology, a real- valued random field X = {X(t), t € M d } is called operator-scaling 
if there exist a linear operator E on M d with positive real parts of the eigenvalues 
and some constant f3 > such that for all constant r > 0, 

{X(r E t), teR d } = {r? X(t), teR d }. (1.2) 

In the above and in the sequel, r E is the linear operator on R d defined by r E = 
Ylri=o ^ nT "rc! E ■ A typical example of Gaussian random fields satisfying (II. 2p is frac- 
tional Brownian sheets introduced by Kamont [13] and other examples have been 
constructed in [U EI]. We mention that (11.20 leads to anisotropy in the "time"- 
variable t, which is a distinct property from those of one-parameter processes. Sev- 
eral authors have proposed to apply such random fields for modeling phenomena in 
spatial statistics, stochastic hydrology and imaging processing (see [5j [21 [8]). 

In this paper, we further extend the notions of operator-self-similarity and operator- 
scaling to multivariate random fields by combining the aforementioned two ap- 
proaches. That is, we will allow scaling of the random field in both "time"-domain 
and state space by linear operators. This is mainly motivated by the increasing 
interest in multivariate random field models in spatial statistics as well as in applied 
areas such as environmental, agricultural, and ecological sciences, where multivariate 
measurements are performed routinely. See Wackernagel [30], Chiles and Delfiner 
[6] and their combined references for further information. We also believe that the 
random field models constructed by Zhang [34] , Gneiting, Kleiber and Schlather 
[11] . Apanasovich and Genton pQ are locally operator-self-similar and their tangent 
fields are operator-self-similar in the sense of Definition 11.11 below. This problem 
will be investigated in a subsequent paper. 
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Throughout this paper, let X = {X(t),t G R d } be a random field with values in 
R m , where d > 2 and m > 2 are fixed integers. In the probability literature, R d is 
often referred to as the "time" -domain (or parameter space), R m as the state space 
and X a (d, m)-random field. We will be careful not to confuse the terminology with 
the space-time random fields in geostatistics. 

The following definition is a natural extension of the wide-sense operator-self- 
similarity and operator-self-similarity in Sato |28| for one-parameter processes to 
(d, m)-random fields. 

Definition 1.1 Let E be a d x d matrix whose eigenvalues have positive real parts. 
A (d,m) -random field X = {X(t),t G R d } is called wide-sense operator- self- similar 
(w.o.s.s.) with time-variable scaling exponent E, if for any constant r > there 
exist an m x m matrix B{r) (which is called a state space scaling operator) and a 
function a r {-) : R d -> R m (both B(r) and a r (-) are non-random) such that 

{X(r E t), t£R d } = {B(r)X(t) + a r (t),t G R d }. (1.3) 

If, in addition, a r (t) = 0, then X is called operator- self- similar (o.s.s.) with scaling 
exponent E. 

Remark 1.1 Here are some remarks about Definition 11.11 

(i) If a random field X is w.o.s.s., then the consistency in (jl.3p implies 

B(nr 2 ) = B( ri )B(r 2 ) = B(r 2 )B(n), Vri,r 2 >0 (1.4) 
and for all n, r 2 > and t £ R d , 

a rir2 (t) = B(n)a r2 (t) + a ri (rf t) = B(r 2 )a ri (t) + a r2 (rft). (1.5) 

(ii) One can also define operator-self-similarity for random fields by extending the 
analogous notion in Hudson and Mason [12] . Namely, we say that a (d, m)- 
random field X = {X(t),t G R d } is operator-self-similar (o.s.s.) in the sense of 
Hudson and Mason with time- variable scaling exponent E, if for any constant 
r > there exist an m x m matrix B(r) and a vector a{r) G R rn such that 

{X(r E t), teR d } = {B{r)X(t) + a(r),t G R d }. (1.6) 

Since the function a(r) does not depend on t G R d , (|1.6p is stronger than 
w.o.s.s. in Definition ll.il but is weaker than the operator-self-similarity. 

Recall that a probability measure [i on R m is full if its support is not contained in 
any proper hyperplane in M m . We say that a (d, m)-random field X = {X(t), t G R d } 
is proper if for each t ^ 0, the distribution of X(t) is full. Then one can verify (see 
e -g- (121 p. 282]) that for a proper w.o.s.s. random field, its space-scaling operator 
B(r) must be nonsingular for all r > 0. 
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We remark that proper w.o.s.s. random fields are special cases of group self- 
similar processes introduced by Kolodyriski and Rosiiiski |14j and can be studied 
by using their general framework. To recall their definition, let G be a group of 
transformations of a set T and, for each (g, t) G G x T, let C(g, t) : R m — > M m be a 
bijection such that 

C(si<72,*) = C{ gi ,g 2 (t))oC(g2,t), \/ gi ,g 2 G G and t G T, 

and C(e,t) = I. Here e is the unit element of G and / is the identity operator 
on M. m . In other words, C is a cocycle for the group action (g,t) *— > g{t) of G on 
T. According to Kolodyhski and Rosihski [14] , a stochastic process {X(t),t G T} 
taking values in M m is called G -self- similar with cocycle C if 

{X(g(t)) ,t€T}± {C(g, t)X(t), t G T}. (1.7) 

Now we take T = R d and G = {r E : r > 0} which is a subgroup of invertible linear 
operators on M. d . It is clear that if a proper (d, m)-random field X = {X(t),t G 
R rf } is w.o.s.s. in the sense of Definition II. 1\ then it is G-self-similar with cocycle 
G, where for each g = r E G G and t G M d , C(g,t) : M. m -> R m is defined by 
C(g,t)(w) = B{r)w + a r (t). Note that C(g,t) is a bijection since X is proper; and 
it is a cocycle because of (ll.4p and (I1.5P . 

In [14] . Kolodyhski and Rosihski consider a strictly stable process X = {X(t), t G 
T} with values in M m which is G-self-similar with cocycle C and characterize the 
minimal spectral representation of X (which is a kind of stochastic integral repre- 
sentation and always exists for strictly stable processes) in terms of a nonsingular 
action L of G on a measure space (S 1 , B(S),fj,), where S is a Borel subset of a Polish 
space equipped with its Borel <r-algebra B(S) and /i is a u-finite measure, and a 
cocycle c : G x S — > {— 1,1} relative to L (see Section 3 of [H] for details). They 
also construct strictly stable processes which are G-self-similar with cocycle G by 
using nonsingular actions L of G on S and {—1, l}-valued cocycle c relative to L (see 
Section 4 of |14j). Their general framework provides a unified treatment for stochas- 
tic processes with various invariance properties (such as stationarity, isotropy, and 
self-similarity) and is particularly powerful when combined with methods from er- 
godic theory to study probabilistic and statistical properties of G-self-similar strictly 
stable processes. See, Rosihski [22], [23] , Roy and Samorodnitsky [23] and Samorod- 
nitsky [26j for recent results on stationary stable processes and random fields. It 
would be very interesting to pursue further this line of research for o.s.s. or more 
general G-self-similar stable random fields. 

The main objective of the present paper is to characterize the permissible forms 
for the state space scaling operator (or simply the space-scaling operator) B(r), 
which provides corresponding information on the cocycle C(g,t). We will also con- 
struct two types of proper o.s.s. symmetric a-stable (d, m)-random fields by using 
stochastic integrals of matrix- valued deterministic functions with respect to vector- 
valued symmetric a-stable (SaS) random measures. Our construction method is 
somewhat different and less general than that of Kolodyhski and Rosihski [14] who 
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use stochastic integrals of real- valued deterministic functions with respect to a real- 
valued strictly stable random measure and who only require their deterministic in- 
tegrands to satisfy certain recurrence equation involving a non-singular action L 
of G on S and a cocycle c : G x S — > {— 1,1} relative to L. See Proposition 4.1 
in [14] for details. The deterministic integrands in our constructions are given in 
terms of O-homogeneous functions (see Definition 2.6 in [3] or Section 2 below). 
Hence the resulting o.s.s. stable (d, m)-random fields in this paper are natural mul- 
tivariate extensions of the familiar linear and harmonizable fractional stable fields. 
To explore the connections between these o.s.s. stable random fields and the G- 
self-similar stable random fields in Proposition 4.1 of Kolodyhski and Rosihski [14] . 
we determine the non-singular action L of G = {r E ,r > 0} on the measure space 
(R d , B(R d ), A d ) and the cocycle c:GxR d ^ {-1, 1} relative to L for the o.s.s. SaS 
random fields constructed in Theorems 12.51 and 12.61 These preliminary results may 
be helpful for applying the powerful tools developed in Rosihski [22j [23] to study 
operator-self-similar SaS random fields. 

The rest of this paper is divided into three sections. In Section 2, we provide 
some preliminaries and state the main results of this paper. Theorem 12.11 proves 
that, under some standard conditions, the space-scaling operator B(r) in fjl .3|) must 
be of the form B(r) = r D for some D G M(R m ), which will be called the state 
space scaling exponent (or the space-scaling exponent) of X. Theorem 12.21 is an 
analogous result for o.s.s. random fields. Theorems 12 . 51 and 12.61 provide general ways 
for constructing proper moving-average-type and harmonizable-type o.s.s. stable 
(d, m)-random fields with prescribed operator-self-similarity exponents. We also 
describe the connection between these random fields and the G-self-similar stable 
random fields in [T4]. In Section 3 we characterize the forms of the space-scaling 
operators and prove Theorems 12.11 and 12.21 The proofs of Theorem 12.51 and 12.61 are 
given in Section 4. It will be clear that the arguments in Hudson and Mason [12] , 
Maejima and Mason [IT] and Bierme, Meerschaert and Scheffler [4] play important 
roles throughout this paper. 

We end this section with some notation. For any integer n > 1, we use A n to 
denote the Lebesgue measure on R n and £>(R™) the Borel algebra. The Euclidean 
norm and inner product in R n are denoted by | • | and (•,•), respectively. Let End(R n ) 
be the set of all linear operators on R n or, equivalently, n x n matrices. The set of 
invertible linear operators in End(]R n ) is denoted by Aut(R n ). Let QiW 1 ) be the set 
of A € Aut(R n ) such that all eigenvalues of A have positive real parts. Let M(R n ) 
be the set of A £ End(R n ) such that all eigenvalues of A have nonnegative real parts 
and every eigenvalue of A with real part equal to zero (if it exists) is a simple root 
of the minimal polynomial of A. 

We will use Co, C\, C2, ■ ■ ■ to denote unspecified positive finite constants which 
may not necessarily be the same in each occurrence. 

2. Main results 

Throughout this paper, E S Q(Mr) is a fixed d x d matrix. E* is the adjoint of 
E; and a € (0, 2] is a constant. 
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Our first result characterizes the form of the space-scaling operator B{r) for a 
w.o.s.s. random field. 

Theorem 2.1 Let X = {X(t),t G R d } be a stochastically continuous and proper 
w.o.s.s. random field with values in R m and time-variable scaling exponent E G 
Q(R d ). There exist a matrix D G M(R m ) and a function b r (t) : (0, oo) x R d -> R m 
which is continuous at every (r,t) G (0, oo) x R d such that for all constants r > 

{X(r E t), teR d } = {r D X(t) + b r (t), teR d }. (2.1) 

Furthermore, X(0) = a a.s. for some constant vector a G R m if and only if D G 
Q(R m ). In this latter case, we define bo(t) = a for all t G R d , then the function 
(r,t) i — y b r (t) is continuous on [0, oo) x R d . 

The operator D will be called the state space scaling exponent (or space-scaling 
exponent). For a given time- variable scaling exponent E G Q(R d ), the correspond- 
ing exponent D may not be unique. In order to emphasize the roles of the linear 
operators E and D, we call X w.o.s.s. with exponents (E,D), or simply [E,D)~ 
w.o.s.s. By combining Theorem 12.11 with Lemma 2.4 in [28], we derive readily the 
following corollary. Of course, (|2.2p also follows from (jl.5p . 

Corollary 2.1 Under the conditions of Theorem \2.1\ the function b r (t) is uniquely 
determined by E and D. Furthermore, 

b rir2 (t) = b ri (rf t) + rf b r2 (t) = b r2 (rft) + r 2 D 6 n (t) (2.2) 

for all n, r 2 > and t G R d . 

The following corollary expresses the function b r (t) in terms of a function of t and 
the scaling exponents E and D. 

Corollary 2.2 Under the conditions of Theorem \2.1l there exists a continuous func- 
tion b(-) : R d \{0} —> R m such that the function b r (t) satisfies 

b r (t) = b(r E t) - r D b{t), Vr > and t G R d \{0}. (2.3) 

If, in addition, D G Q(R m ) and X(0) = a a.s., where a G R m is a constant vector, 
then we can extend the definition of b(-) to R d by defining 6(0) = a such that \2.3\) 
holds for all r > and t eR d . 

The proof of this corollary is based on the polar coordinate representation of 
t G M' i \{0} under operator E given in [H p. 317] (the definition is recalled below) 
and will be given in Section 3. 

The next theorem is an analogue of Theorem 12.11 for o.s.s. random fields. 

Theorem 2.2 Let X = {X(t),t G R d } be a stochastically continuous and proper 
random field with values in R m . 
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(i) If X is o.s.s. with time-variable scaling exponent E G Q(R ), then there exists 
a matrix D G M(R m ) such that for all r > 

{X(r E t), teR d } = {r D X(t), t G R d ). (2.4) 

Moreover, D G Q(P m ) i/ and onZy i/Jf(0) = a.s. 

(it) If X is o.s.s. with time-variable scaling exponent E G Q(R. d ) in the sense of 
Hudson and Mason, then there exist a matrix D G AI(R m ) and a continuous 
function b{r) : (0, oo) — > W n such that for all constants r > 

{X(r E t), t£R d } = {r D X(t) + 6(r), teR d }. (2.5) 

A (d, m)-random field X = {X(t),t G R d } is called operator-self-similar with 
exponents (E,D) (or (E, D)-o.s.s.) if (12. 4p holds. By Corollary 12.21 we see that if 
X is a w. o.s.s. (d, m)-random field as in Theorem 12 J] then the (d, m)-random field 
Y = {X(t) —b(t),t G M rf \{0}} is operator-self-similar with exponents (E,D). Using 
the terminology of Sato [28], we also call the function b(t) in Corollary 12.21 the drift 
function of (d, m)-random field X. 

Recall that a (d, m)-random field X is said to have stationary increments if for all 
h G R d , 

{X(t + h)- X(h), t G R d } = {X(t), teR d }. (2.6) 

Now we turn to construction of interesting examples of stable o.s.s. (d, m)-random 
fields with stationary increments, by using stochastic integrals with respect to a 
stable random measure. We refer to Samorodnitsky and Taqqu [27] for a systematic 
account on the latter. For simplicity we will only consider symmetric a-stable (SaS) 
random fields and the main idea comes from [4], [19] and [17]. By using stochas- 
tic integral with respect to a strictly stable random measure one can extend the 
construction to obtain strictly stable o.s.s. (d, m)-random fields. Kolodyhski and 
Rosihski [J3] use this more general approach. 

For any given operators E G Q(R d ) and D G Q(R m ) we construct (E, D)-o.s.s. 
a-stable random fields by using stochastic integrals with respect to a symmetric 
a-stable random vector measure (when a = 2 the resulting o.s.s. random fields are 
Gaussian). For this purpose, we recall briefly the definitions of stochastic integrals 
with respect to vector- valued a-stable random measures. 

Let (fi, P) be the underlying probability space and let L°(Q) be the set of all 
R m - valued random vectors defined on (f2, P). Let S m -i be the unit sphere in R m 
with the Borel algebra B(S m -i). 

Let K be a a-finite measure on R d x S m -\ such that for any A G B(R d ), K(A x •) 
is a symmetric finite measure on (5 m _i, B(S m ^i)). Denote 

M := {A G B{R d ) : K{A,S m - X ) < oo}. 

We first give the definition of a vector-valued symmetric a-stable (SaS) random 
measure. 
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Definition 2.1 An W 11 -valued SaS random measure on (R d ,B{R d )) with control 
measure K is an independently scattered a-additive M. m -valued set function M : 
M — > L°(tl) such that, for every A G M, the random vector (M 1 (A), • • • , M m (A)) 
is jointly SaS with spectral measure K(A,-). Here, the meaning of "independently 
scattered" and "a-additive" is the same as in Section 3.3 of [Wi 



One can apply Kolmogorov's extension theorem to show that R m -valued SaS 
random measure M in Definition 12.11 exists, with finite-dimensional distributions 
characterized by 



Eexp I {^{e^MiAj)) 



{/■ r- km 



K{dx,ds)}, (2.7) 



where Aj G M and Oj = (0 jt i, ■■■ , 9 j>m ) G R m for all k > 1 and j = 1, • • ■ , k. 

In this paper, unless stated otherwise, the control measure K will always be as- 
sumed to have the form K(A, B) = X d {A)T(B) for all A G B(R d ) and B G B(S m -i), 
where is the Lebesgue measure on W d and T(-) is the normalized uniform measure 
on 5 m _i such that for all 6 = (0i, • • • , m ) G M m , 



m 

i=i 



r(ds) = \e\ 



Therefore, for disjoint sets Aj G Ai, j 

k 



1, 2, • • • , k, Eq. ()2.7|) can be written as 



Eexp U^2,{6j, M{Aj))\ = exp 

7 = 1 ' 



k 

E 



(2.8 



For any real mxm matrix Q, let ||Q|| := max| x .| =1 \Qx\ be the operator norm of Q. 
It is easy to see that for Q u Q 2 G End(M m ), ||QiQ 2 || < IIQill ■ ||Q2||- The following 
theorem is an extension of Theorem 4.1 in [T7] and defines stochastic integrals of 
matrix-valued functions with respect to a vector-valued SaS random measure. 

Theorem 2.3 Let {Q(u),u G M. d } be a family of real m x m-matrices. If Q(u) is 
B(M. d ) -measurable and J„ d \\Q(u)\\ a du < oo, then the stochastic integral 



i(Q) ■■- 



Q(u)M(du) 



is well defined and it is a symmetric a-stable vector in M m with characteristic func- 
tion 



E 



exp 



\Q( u ye\ a du^, 



ye g 



(2.9) 
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It follows from (|2.9p and Lemma 13.21 below that if the matrix Q(u) is invertible 
for u in a set of positive A^-measure, then the distribution of I(Q) is full. This fact 
is useful for constructing proper SaS random fields. 

One can also define stochastic integrals of complex matrix-valued functions with 
respect to a complex vector-valued SaS random measure M defined as follows. Let 
M be an M 2m - valued SaS-random measure on {R d ,B(R d )) with control measure K = 
\d x r, where V is the normalized uniform measure on S^m-i- Define the complex- 
valued ^SaS-random measures Mk = + iM m+ k for all k = 1, • • • , m. Then 
M = (M\, • • • , M m ) is a C m -valued SaS-random measure with control measure K. 
Its real and imaginary parts are Mr = (M 1; • • • , M m ) and Mj = (M m+1 , • • • , M 2m ), 
respectively. The following theorem defines stochastic integrals of complex matrix- 
valued functions with respect to M. 

Theorem 2.4 Let {Q±(u), u G W 1 } and {(^(it), u € M. d } be two families of real 
m x m-matrices. Let Q(u) = Qi(u) + iQ2(u) for all u G M. d . If Q\(u) and Q2(u) 
are B(M. d ) -measurable and f^d(\\Qi(u)\\ a + ||Q2(ii)|| a ) du < oo, then 



I(Q) := Re [ Q(u) M(du) 



is well defined and it is a symmetric a-stable vector in M m with its characteristic 
function given by 



E 



exp i - J ^ (J \Qi(u)*9\ 2 + \Q 2 (u)*9\ 2 \ du\, V9 G 



(2.10) 



It follows from fl2. lQf) and Lemma 13.21 below that if the matrix Q\(u) or Q2{u) is 
invertible for it in a set of positive A^-measure, then the distribution of I(Q) is full. 

Based on the above stochastic integrals, we can construct moving-average type 
or harmonizable-type a-stable random fields by choosing suitable functions Q and 
Q. In order to obtain o.s.s. random fields, we will make use of the 0-homogeneous 
functions and the (/3, 0)-admissible functions as in [1]. 

Suppose G Q(U. d ) with real parts of the eigenvalues < a\ < a 2 < ■ ■ ■ < a p 
for p < d. Let q denote the trace of ©. It follows from [3J p. 314] that every 
x G M. d \ {0} can be written uniquely as x = t(x) & 1(x) for some radial part t(x) > 
and some direction l(x) G So such that the functions x i-> t(x) and x i— >• l(x) are 
continuous, where So = {x G t(x) = 1}. It is well-known that t(x) = t(—x) 
and r(r e x) = rr(x) for all r > and x G M. d . Moreover, So is compact; t(x) — > oo 
as \x\ — > oo and t(x) — > as |x| — > 0. In addition, Lemma 2.2 in [3] shows that 
there exists a constant Cq > 1 such that for all x, y G M d 

T(x + y) <C (T(x) + T(y)). (2.11) 

For convenience, we call (t(x),1(x)) the polar coordinates of x under operator 0. 
According to Definition 2.6 in [4], a function (j) : M. d — > C is said to be ©-homogeneous 
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if 0(r e x) = rcj)(x) for all r > and x € \ {0}. Obviously, if <p is ©-homogeneous, 
continuous on M. d and takes positive values on R d \ {0}, then 0(0) = 0, 

= max (f)(9) > and = min (f)(9) > 0. (2.12) 

Let /3 > 0. Recall from Definition 2.7 in [1] that a function t/> : R d — > [0, oo) is 
called (/3, 0)-admissible, if ^(s) > for all x / and for any < A < B there 
exists a positive constant C\ > such that, for A < \y\ < B, 

t(x)<1 \^(x + y)-^(y)\<C l T(xf. 

For any given matrices E £ (5(R d ) and D € (^(IR™ - ), Theorem 12.51 provides a class 
of moving-average-type o.s.s. a-stable random fields with prescribed self-similarity 
exponents (E,D). 

Theorem 2.5 Suppose (f> : M. d i— > [0,oo) is an E -homogeneous, ((3, E)- admissible 
function for some constant (3 > 0. Let q be the trace of E, H be the maximum of 
the real parts of the eigenvalues of D € Q(W n ) and let I be the identity operator in 
W 71 . If H < f3, then the random field 

X^(x) = [ \(f)(x - y) '" 1 / - 0(- y ) B -^/«j M(dy), x&R d (2.13) 



is well defined, where the stochastic integral in A2.13\) is defined as in Theorem W. 
Furthermore, X§ = {X ( / ) (x),x 6 R d } is a stochastically continuous (E, D)-o.s.s. 
SaS-random field with stationary increments. 

Remark 2.1 We can choose E and D to ensure that the SaS-random field X is 
proper. A sufficient condition is that q/a is not an eigenvalue of D. This implies 
that, for every x € R d , the operator <f>(x — y) D ~i I / a — (f>(—y) D ~ qI / a is invertible 
for y in a subset of R d with positive Lebesgue measure, which ensures that the 
distribution of X ( / ) (x) is full. 

When m = 1 and D = HI, Theorem 12.51 reduces to Theorem 3.1 in Bierme, 
Meerschaert and Scheffler [1]. For a general D £ Q(R m ), the following example of 
X^ is instructive. Let E = (ey) be the diagonal matrix in Q(R d ) with ejj = T^ 1 , 
where jj € (0, 1) (1 < j < d) are constants. It can be verified that there exists a 
constant C2 > 1 such that the corresponding radial part t(x) satisfies 

d d 

C 2 X E w 7j - T ^ - ° 2 E w 7j ( 2 - 14 ) 

3=1 j=l 

for all x € R d . Note that the function 4>(x) = Ylj=i \ x jV j ls -^-homogeneous and 
(f3, inadmissible with (3 = 1. This latter assertion follows from (|2.14p and the 
elementary inequality \x + y| 7 < |x| 7 + |y| 7 if 7 € (0, 1). Let D € Q(R m ) be as in 
Theorem [231 then X^ = {X^x), x € R d } defined by 



X^(x) 



d s D-ql/a / d ^ D-ql/a 

j= i / \j=i 



M(dy) 
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is an (E, D)-o.s.s. SaS random field with stationary increments. Moreover, since 
H < 1 and q/a > 1 (we have assumed d > 2 in this paper), we see that is proper. 

Similarly to Theorem 12.51 we can construct harmonizable-type o.s.s. SaS stable 
random fields as follows. 

Theorem 2.6 Suppose ip : M. d i— > [0, oo) is a continuous, E* -homogeneous function 
such that ip(x) ^ for i / 0. Let q be the trace of E and let I be the identity 
operator in M. m . If D G Q(W n ) and its maximal real part of the eigenvalues H < a\, 
where a\ is the minimal real part of the eigenvalues of E, then the random field 

X^x) = Re f (e i{x ' y) - l) ^{y)- D ~ qI/a M(dy), Vx G R d , (2.15) 



is well defined, where the stochastic integral in $2.15\) is defined as in Theorem \2. 
Furthermore, X^ = {X^,{x),x G M. d } is a stochastically continuous, proper (E,D)- 
o.s.s. SaS -random field with stationary increments. 

Remark 2.2 Unlike in Theorem 12.51 Xj, in Theorem 12.61 is always proper. 

Theorem 12.61 is a multivariate extension of Theorem 4.1 and Corollary 4.2 of 
Bierme, Meerschaert and Schemer [4]. To give a representative of the harmonizable- 
type o.s.s. in Theorem 12.61 again we take E = (e^) G Q(M d ) to be the diagonal 
matrix as above. Let ip(x) = Y^j=i l ;E jl 7j > which is E* -homogeneous. Then, for 
any D G Q(R m ) with its maximal real parts of the eigenvalues H < min{7~ 1 }, the 
SaS-random field X^ = {X 1 p(x), x G M. d } defined by 

^ {x) ' Re L (T ^Z-) D ^ mv) (2,16) 

is proper and (E, D)-o.s.s. with stationary increments. In the special case of D = I, 
the stable random field X^ has been studied in Xiao [32]. We believe that the 
argument in proving Theorem 3.4 in [32] can be applied to show that X^ has the 
property of strong local nondeterminism, which is useful for establishing the joint 
continuity of the local times of X^. 

The o.s.s. SaS (d, m)-random fields in Theorems I2.5I and I2.6I provide concrete 
examples for the G-self-similar stable random fields in Proposition 4.1 of Kolodyiiski 
and Rosihski [14 . Recall that the o.s.s. SaS random fields in Theorems 12.51 and 
12.61 are G-self-similar with cocycle C, where G = {r E ,r > 0} and C{r,t) = r D 
for every r > and t G M. d . In the following we provide non-singular actions 
of G = {r E ,r > 0} on (R d ,B(R d ), X d ) and cocycles c : G x R d -> {-1,1} (or 
{zGC:|z| = l}in the complex case) such that the integrands in (|2.13p and (|2.15p 
satisfy the recurrence equation (4.1) in Kolodyhski and Rosihski [14]. 

For the o.s.s. SaS random field X^ in Theorem 12.51 the non-singular action of G 
on R d is L r (s) = r E s, and the cocycle c(r, x) = 1. A change of variable shows that 

dArf 
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where q is the trace of E. By using (|2.17p and the .E-homogeneity of 4> one can 
verify that the family of integrands {f x , x G R d } in Theorem 12,51 where 

Uy) = <f>(x - yf-* 1 /* - H-y) D ~ qI/a 

is a matrix-valued function, satisfies 

f rEx {y) = c(r, L r -iy)| d(Ad d ° A ^' l) Y "cir, x)f x o L r -i (y), Vy G M d , (2.18) 
which is an analogue of the recurrence equation (4.1) in Kolodyhski and Rosihski 

For the o.s.s. SaS random field in Theorem 12.61 the non-singular action of 
G on R d is L r (s) = r E * s and the cocycle c(r,x) = 1. Then, by using (I2.17|) and 
the £7*-homogeneity of ifi one can verify that the family of integrands { f x , x G R rf }, 
where 

Uy) = {^ v) -i)m- D - qI/a , 

satisfies the recurrence equation (12.181) with L being replaced by L. 

3. Characterization of space-scaling exponents: Proofs of Theorems 12.11 
and[2H 

In this section, we prove Theorem 12.11 The main idea of our proof is originated 
from pL2j and [28]. We will make use of the following lemmas which are taken from 
[28] and [29], respectively. 



Lemma 3.1 (j28l Lemma 2.6]) 

For any integer n > 1, H G Q(M. n ) if and only if lim^o r H x = for every x G M n . 
H G M(W l ) if and only if limsup r | \r x\ < oo for every x G M n . 

Lemma 3.2 ( [SHI Proposition 1]) 

A probability measure fx on IR n is not full if and only if there exists a vector y G 
M n \{0} such that \jl{cy)\ = 1 for all c£l, where fl is the characteristic function of 
li. 

For r > and E G Q(M d ) fixed, define G r to be the set of A G Aut(M m ) such that 
{X(r E t), t €R d } = {AX(t) + b(t), t £R d }, for some function b : R d -> R m . Let 

G = Ur>0 

Lemma 3.3 The set G is a subgroup of Aut(R m ) . In particular, the identity matrix 
I G G\; A G G r implies A~ l G Gu r ; A G G r and B G G s imply AB G G sr . 

Proof. This can be verified by using the above definition and the proof is elementary. 
We omit the details here. □ 

Lemma 3.4 The following statements are equivalent: 

(1) There exist a sequence {r n ,n > 1} with r n J, and A n G G rn such that A n tends 
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to A G Aut(R m ). 

(2) {X(t), t G R d } = {X(0) + (f>(t), t G R d }, where <ft is unique and continuous on 
R d . 

(3) G = G S for all s > 0. 

(4) G s n G r ^ for some distinct s, r > 0. 

Proof (1)=K2). Assume (1) holds then we have that {X(r E t), teR d } = {A n X(t) + 
b rn (t), t G R d }. By Lemma |3. II and the stochastic continuity of X, we derive that 

there is a function b(t) such that {X(t), t€R d } = {A^X^-A^b^), t G R d } and, 

in particular, X(0) = A^X(0) - A^tyO). This yields (2) with <p(t) = A _1 6(0) - 
A~ l b(t). The continuity of follows from the stochastic continuity of X and the 
uniqueness of </> follows from Lemma 2.4 in [28]. 

(2) =K3) Suppose (2) holds and A G G P . Then 

{X(0) + 0(r B t), t G M d } = {X(r s i), (GK d } = {AY(t) + 6 r (t), t G 
Hence for all positive numbers s ^ r, 

{X{s E t), t G R d } = {X(0) + 0(a B t), t G R d } 

= {AX{t) + b r {t) - cj)(r E t) + 0(s B i)> t G R d }. 

Thus yl G G s , which shows G r C G s . By symmetry, we also have G s C G r . Therefore 
G r = G s for all s ^ r, and hence G r = G. 

(3) =>(4) This is obvious. 

(4) =>(1) Now we assume (4) holds for some s < r. Let A G G s C\ G r . Since 

{X(s £ t), £ G R d } = {AX(t) + b s (t), t G M d } and {X{r E t), t £ R d } = {AX(t) + 

6r(t), t G we obtain that {X(s E t), t<ER d } = {X(r E t) + ip(t), t G R d } for some 
function ip : R d ->■ M m . Then 

{X((s/r) s t) ; t G M d } = {X(i) + ^(r^i), i G (3.1) 
This shows that I G G s / r . Let c n = (s/r) n . By iterating (|3.ip we derive that 

{X(c E t), t G R d } = {X(t) + ^n(t), * G 

n-l 

where ip n (t) = Yl 4>( c f r~ E t). Hence / G G Cn for all n > 0. bmce c n ->■ and I G 

8=0 

Aut(M m ), we arrive at (1). □ 

Lemma 3.5 Assume G ^ G s for some s > 0. If A n G GV n , A G Aut(R m ) and 
A n —¥ A as n — > oo, £/ien i/ie sequence {r n } converges to some r > as n -4 oo and 
A g G r . 

Proof. Suppose that {r nfc } is a subsequence of {r n } and that {r nk } converges to 
some r G [0, oo]. Then < r < oo. In fact if r = 0, then A nk —¥ A and Lemma 
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13,41 imply G = G s for all s > 0, which is a contradiction to the assumption. On the 
other hand, if r — > oo, then A~ k G G r -i — > A -1 and r~ k -4 0. By Lemma 13,41 we 
also get a contradiction. It follows from 

{X(rlt), (GK^ {^ fe X(t) + b rnk (t), t G M d } 
and the stochastic continuity of X that 

{X(r^), i G M. d } = {AY(t) + 6 r (i), t G M d } 

for some function b r . Therefore A £ G r and hence from Lemma 13.41 we infer that 
all convergent subsequences of {r n } have the same limit r. Consequently, {r n } 
converges to r > 0. □ 

From Lemma 13.41 and Lemma |3.5[ we derive the following result. 
Corollary 3.1 If G ^ G s for some s > 0, then Gi is not a neighborhood of I in G. 

Proof. From Lemma 13.41 the assumption that G ^ G s for some s > implies 
G r n G\ = for all r 7^ 1. Therefore, to prove the corollary, it is enough to show 
that there exists a sequence A n £ G Tn such that r n ^ 1 and A n — > / as n — > oo. 
This can be proved as follows. 

Let {r n } be a sequence with r n ^ 1 and r n — >• 1 as n — > oo. Take S n € G r „- Then 
by the convergence of types theorem (see, e.g., |29t p. 55]), {B n } is pre-compact in 
Aut(M m ). Hence we can find a subsequence {B nk } such that B nk — > B € Aut(M m ). 
By Lemma 13.51 we have B £ G\ and thus -B" 1 E Cri- Furthermore, by Lemma 13.31 
G r „ k 3 B B nk — > I € Gi. Let A^ = B~ 1 B 1lk , then the sequence {Ak}k is what we 
need. □ 

Using the above results, we give the proof of Theorem 12.11 as follows. 

Proof of Theorem \2.1\ From Lemma 13.41 we only need to consider two cases. 

Case 1: G = G s for all s > 0. By Part (2) of Lemma [331 we derive that for all 
constant c > 0, 

{X(r E t), t£R d } = {X(0) + cj)(r E t), t G R d } 

= {X(t) + ^{r E t)-^(t),teR d }. 

Hence f|2. 1|) holds with D = 0, which is the matrix with all entries equal 0, and 
b r (t) = <f>(r E t) - 4>{t). 

Case 2: {G s ,s > 0} is a disjoint family. In this case, G is a closed subgroup of 
Aut(M m ). Define ??: G -> E by = Ins if A G G s . It is well-defined and, from 
Lemma 13.31 and Lemma 13. 5( is a continuous homomorphism between the group G 
and the group (R, +). Let T(G) be the tangent space to G at the identity /. It is 
well-known that the image of T(G) under the exponential map is a neighborhood of 
the identity of G; see [TJ p. 110.]. Therefore, by Corollary 13. 1| there exists A G T(G) 
such that e A g" G\. Furthermore, by the same arguments used in the proof of 
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Theorem 2.1 of [El p. 288], we know there is a D G End(M m ) such that G G s for 
every s > 0. This implies that 

{X{r E t), t G M. d } = {r D X(t) + 6 r (i), i G R d }. (3.2) 

for some function b r (t). Note that the linear operators r E and r D are continuous 
on r G (0, oo). By the convergence of types theorem, it is not hard to see that 
b r (t) is continuous in (r, t) G (0, oo) x ~R d . In order to verify the fact D G M(M. m ), 
we let {X Q (t),t G M d } be the symmetrization of {X(t),t G M d } and let /x(t) be the 
distribution of X (t). Then by (f3T2j) 

A*(r B t) = r D M (t), 

for all r > and i G M d . Therefore, the characteristic function of [i(t), denoted by 
fit(z) (z G R m ), satisfies 

jl r E t (z) = fi t (r D * z) (3.3) 

for every r > and t G R d , where D* is the adjoint of D. Suppose D G" M(M. m ), then 
D* <=L M(R m ) either. By Lemma l3Tj we can find r n — > and z G M m such that 
zq\ — > oo. Let a n = zo| _1 - Then by choosing a subsequence if necessary, we 
have that a n r®* zq converges to some z\ G M m with \z\\ = 1. From ()3.3j) . it follows 
that for all c G M 

J2 r E t (ca n z ) =fit(ca n r^* z ). (3.4) 

Letting n — > oo, since Lemma 13.11 implies r E t — > 0, by the continuity of £tt(-), we 
have that Jlt(czi) = /Uo(0) = 1 for all cel. It follows from Lemma 13.21 that X(t) 
is not full in ]R m . This contradicts the hypothesis that X is proper. Consequently, 
the matrix D in (|3.2I) belongs to M(M. m ) and the function b r (t) is continuous in 
(0, oo) x R d . 

Now we prove that X(0) = a a.s. for some constant vector a G M m ) if and only 
if D G Q(M m ). From Lemma 13. 1\ it can be shown that, if X is a stochastically 
continuous w.o.s.s. random field and D G Q(R m ), then ^(0) = const, a.s. Con- 
sidering the converse assertion, we note that, in this case, the symmetrization of 
{X(t),t G R d }, i.e. {X (t),t G R d }, satisfies X (0) = a.s. If D & Q(M. m ), then by 
Lemma [3.11 we can find r n — > and zq such that |r^*2;o| does not converge to 0. Let 
a n = \r E zq\~ 1 . Then choosing a subsequence if necessary, by the fact D G M(M m ), 
we have that a n converges to a finite a > and that a n r^ zq converges to some 
z\ G M m with \zx | = 1. By using (|2.ip and the same argument as that leads to (|3.3[) 
and (13. 4p we derive 

ji r E t (ca n z ) = Jl t (ca n r E ' z ) (3.5) 

for all c G M. Letting n — > oo, we have that /xj(c^i) = /^(co^o) = 1- Then by Lemma 
[331 X(t) is not full in R m . This contradiction implies D G Q(M m ). 

The last assertion follows from the stochastic continuity of X and (|2. 1 j) . This 
finishes the proof of Theorem 12.11 □ 



15 



16 



Proof of Corollary \2.2\ For every t G R d \{0} we use polar coordinate decomposition 
under the operator E to write it as t = TE{t) E l{i). We define b(t) = b TE ^(l(t)^ for 
t G R d \{0}. Then from {2J2]) we derive that for all r > and t G R d \{0}, 

b rTE{t) {l(t)) = b r (r E (t) E l(t)) +r D b TE{t) (l(t)), 

which can be rewritten as 

b{(rr E (t)) E l(t))=b r {t)+r D b{t). 

This implies b r (t) = b(r E t) - r D b{t) for all r > and t G R d \{0}. In the case when 
X(0) = a a.s., (|2.ip implies b r (0) = a — r D a, which shows that f|2.3j) still holds for 
t = 0. □ 

Proof of Theorem \2.2\ The proof is similar to the proof of Theorem 12. 1\ with some 
minor modifications. For proving Part (i), we define G r to be the set of A G Aut(R m ) 

such that {X(r E t), t £ R d } = {AX(t), t £ R d }; and for proving Part (ii), we define 

G r to be the set of A G Aut(M m ) such that {X(r E t), t G R d } = {AX{t) + 6(r), t G 
for some function b : (0, oo) — > W 71 . The rest of the proof follow similar lines 
as in the proof of Theorem 12.11 and is omitted. □ 

We end this section with two more propositions. Proposition 13.11 shows that, if a 
(d, m)-random field X is w.o.s.s. with time-variable scaling exponent E, then along 
each direction of the eigenvectors of E, X is an ordinary one-parameter operator- 
self-similar process as defined by Sato [28J. Proposition ^. 21 discusses the relationship 
between w.o.s.s. random fields and o.s.s. random fields in the sense of Hudson and 
Mason (see (ii) in Remark 1.1). 

Proposition 3.1 Let X = {X(t),t G R d } be a stochastically continuous and proper 
(E, D)-w. o.s.s. random field with values in M m . Let X be a positive eigenvalue of 
E and £ G R d satisfy E£ = A£. Denote b r (u) = b r (u^) for all u G R. Then the 
following statements hold. 

(i) There exists a continuous function f{u) from R\{0} to R m , such that b r (u) = 
f(ur x ) — r D f(u) for all u ^ and r > 0. 

(ii) If D G Q(R m ), then f(u) can be defined at u = such that f(u) is continuous 
in R. Moreover, the stochastic process Y = {Y(u), u G R} defined by Y(u) = 
X(u£) — f(u) satisfies that for any r > 

{Y(ru), u G R} = {r D/x Y(u), u G R}. 

Proof. By Corollary I2.H we have that 

iVir 2 «) = b ri {r E u£) + r E b r2 (uC) 
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for all n, T2 > 0. Since £7£ = A£ and r 2 u^ = ur A £, we have 

b rir2 «) = b ri (r A <) + rf 6 r2 «) . (3.6) 

Define /(ti) = & m i/a(£) for u > and /(ti) = &i u ii/a ( — £) f° r u < 0. Then the 

continuity of f{u) on R\{0} follows from the continuity of b r (t). Moreover, from 
(ESI) it follows that 



_ brM) = -rf/(r 2 A ) + /(r A r 2 A ), (3-7) 
&ri H 2 A ) = -*f /Hi) + /HM)- (3-8) 
Writing u = r 2 or — r 2 and r = n, we see that (13.7|) and fj3.8f) yield that 

& r («) = /(ur A )-r D /(u) (3.9) 

for all r > 0, u 7^ 0. This proves (i). 

Suppose D G Q(R m ). Lemma O implies that r D X(£) -> and r D X(-£) -> in 
probability as r — )• 0. Theorem 12.11 and the convergence of types theorem indicate 
that, as r — > 0+, the limits of b r (£) and b r {— £) exist and coincide. Hence, we can 
define /(0) := lim r _j.o &r (0 • Then f(u) is continuous in R. Combining (|2.ip and 
(l3~9l) yields that for all r > 0, u G R, 

(X(r A <), «£!} = {X(r £ <), iieK} 

= (r D X«) + /(ur A ) - r D /(u), u G R}. 

Hence for the process y = {y(u),ii G R} defined by y(u) = X(u£) — f(u), we have 

{Y(r x u),u G R} = {r D Y(u),u G R}. Equivalently, y is D/A-o.s.s. This finishes the 
proof. □ 

Proposition 3.2 Lei X = {X(t),t G R rf } be a stochastically continuous and proper 
(E, D)-w.o.s.s. random field with values in R' m . Suppose E has two different positive 
eigenvalues \\ and A2. Then X is o.s.s. in the sense of Hudson and Mason if and 
only if b T (t) in \2.V\) only depends on r and \t\ for all r > and t G R rf . 

Proof. The "necessity " part is obvious, because, for every (E, D)-o.s.s. random 
field in the sense of Hudson and Mason, the function b(r) does not depend on t. 

In the following, we prove the sufficiency. Suppose b r {t) only depends on r and 
\t\ for all r > and t G M. d . Then we can find a function g on R 2 such that 
b r (t) = 9( r , \t\)- By Corollary 12 .1\ we have that for all r\,r2 > and t G IR d 

g( n r 2 , \t\) = g(n, \ift\) + rf <?(r 2 , \t\). (3.10) 

Let £1, £2 be the eigenvectors of E corresponding to Ai and A 2 , respectively. Without 
loss of generality, we assume |£i| = I62I = 1 and A2 < Ai. Then from (I3.10p . we have 
that 

g(nr 2 , 1) = g(n, r 2 x ) + rf g(r 2 , 1), 
g(nr 2 , 1) = #(ri,r 2 2 ) + rf g(r 2 , 1), 
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where we have used the facts r E £\ = r Al £i and r E & = f X2 ^,2- Therefore, we derive 
that g(r, u Xl ) = g(r, u^ 2 )for any r > and u > and hence, for all n > 1, 

g(r,u) =g(r,u x ^) = g(r,u X ^ X "). (3.11) 

Note that by Theorem 12.11 g( r i u ) is continuous on (0, oo) x [0, oo). Therefore, 

g(r, 0) = lim g(r, u) (3-12) 

and for any u > 0, letting n — > oo, from (|3.1ip we get that 

g(r,u) =g(r,l). (3.13) 

Combining f)3. 12j) with (|3.13j) . we obtain that g(r, 0) = g(r, 1) and hence for all r > 
and u > 0, g(r,u) = g(r, 1). This means b r (t) = g(r, 1) is independent of t. Hence 
the random field X is o.s.s. in the sense of Hudson and Mason. □ 



4. Construction of o.s.s. stable random fields: Proofs of Theorems I2T3H2. 61 

This section is concerned with constructing (E, D)-o.s.s. random fields by using 
stochastic integrals with respect to SaS random measures. In particular, we prove 
the remaining theorems in Section 2. 

Note that Theorem 12.31 is a multiparameter extension of Theorem 4.1 in [17] and 
can be proved by using essentially the same argument with some modifications. 
Hence the proof of Theorem 12.31 is omitted here. In the following, we first prove 
Theorem E5 

Proof of Theorem \2.4\ We divide the proof into two steps. 
(1) When Q(u) is a simple function of the form 

k k 

Q{u) = Qi( u ) + iQ 2 {u) =^Rjl Aj {u) + i^2ljl Aj (u), (4.1) 

J=l 3=1 

where Rj, Ij € End(M m ) and Aj,j = 1, 2, • • • ,k are pairwise disjoint sets in A4, we 
define 

k 

HQ) = ( R j M R( A i) ~ IjMi(Ajj). 

3=1 

Then for any 9 G W 71 , from (|2.8p . we obtain that 

( k 



E 



e i(0J(Q)) 



exp\ -^(|^| 2 + |/;0| 2 ) Q/2 A(^.) 



3=1 



:expi - J d (\Qi(u)*9\ 2 + \Q 2 (u)*e\ 2 y /2 du\. (4.2) 
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(2). When {Q(u)} fulfills f Rd {\\Qi(u)\\ a + \\Q 2 (u)\\ a )du < oo, we can choose a 

sequence of simple functions {Q^(u) = Qi (u) + iQ^\u)} of the form (|4.1h such 
that as n — > oo, 

/ \\Qi(u)* -Q^iuTfdu^O (4.3) 

and 

>(n) 



||Q 2 («)*-Q^(u)*fdn^0. (4.4) 
By the linearity of /(■) we have 

i(Q [n) ) - HQ (e) ) = W {n) - Q (e) ), 

andE(e i < 9 ' 7 W (n) -« ( ' ) ))) equals 

exp { ~ L (^ ){uT ~ ^M*)*! 2 + I &^ " Q?V)*)<) a/2 du j 

>exp(-/ |(Q^(«)*-QW(u)*)e| d«- / |(^ n) (n)*-Qf( U )*)0| a dnl 



which converges to 1 as I, n -^oo^by P~3j) and (jP|) . Thus i"(Q( n )) - I(Q^) ->■ in 
probability as n — s> oo, and /(QW) converges to an M m -valued random vector in 
probability. It is easy to see that the limit does not depend on the choice of {Q^}- 
Therefore, we can define I(Q) as the limit of I(Q^), and hence 

E (V<MQ)>) = Jim E{e l ^ (n) ^ 

= exp| - J ^ (\J\Qi(u)*9\ 2 + |Q 2 (n)*0| 2 ) Q dn|. 

The proof of Theorem 12.41 is completed. □ 

In order to prove Theorem 12.51 and Theorem 12.61 we will use the following change 
of variable formula from [3]. 

Lemma 4.1 ([4, Proposition 2.3]) Let E € Q(R d ) be fixed and let (t(x), 1{x)) be 
the polar coordinates of x under the operator E. Denote So := {t(x) = 1}. Then 
there exists a unique finite Radon measure a on So such that for all / € L 1 (R d , dx), 

r poo r 

/ f(x)dx= / / f(r E 6)a(de)r q - 1 dr. 

JR d JO JY. 

We also need the following lemma which is due to Maejima and Mason |17| . For 
more precise estimates on ||r D || see Mason and Xiao j!9j . 
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Lemma 4.2 Let D 6 Q(M m ) and let h > and H > be the minimal and maximal 
real parts of the eigenvalues of D, respectively. Then for any 5 > 0, there exist 
positive constants C3 and C4 such that 



C 3 r h ~ s , i/0<r<l, 
C 4 r H+s , ifr>\. 



\r D \\ < 



Now we are in position to prove Theorem 12.51 
Proof of Theorem \2.5l We divide the proof into four parts. 

(i). First we show that the stochastic integral in (|2.13j) is well defined. By Theorem 
12.31 it suffices to show that for all x € W 1 

T°L(x) = I <f>(x - y) D -^/ a - cb(- y ) D -^/- "dy < 00. (4.5) 
Jm. d 

Let (t(x), l(x)) be the polar coordinates of x under operator E. By the fact that <f> 
is .^-homogeneous, we see that 

4>{y) = r(yW(y)) v y e R d . 

Then by (|2.12p . we have that 

m 4 r(y) < <f>(y) < M^r(y). (4.6) 
Therefore, there exists a constant C5 > such that 

||0(l/) £, ~ ,I/a |r<C f 5 ||T(y) IJ -« J / a || a . 

Note that 

Mi = sup ||r fi || > and M 2 = sup \\r E \\ > 

are finite because r E is continuous in r and \\r E \\ 7^ for all r > 0, and that 
< m = inf \y\ < M = sup \y\ < 00, 

since So is compact and So- Therefore, from 

r E (y)y = ^ E (y)r(y) E l(y) = (<p-\y)r(y)) E l{y) 
and (I4.6p . it follows that 

< < \^ E {y)y\ < MM 2 < 00. (4.7) 
1V1\ 

Since <j) is (/?, i?)-admissible, for any z with < \z\ < MM 2 there exists a positive 
constant Ci > such that 

\<f)(x + z) - (j)(z)\ < dT(xf (4.8) 
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for all x € M rf with t(x) < 1. For any 7 > 0, on the set {y € M d : r(y) < 7}, we 



have 



(x - y) 



D-ql/a 



<K-v) 



D-ql/a 



< 2 



(x - y) 



D-ql/a 



+ 2 



<K-v) 



D-ql/a 



Consequently, by Lemma 14. 1\ Lemma 14.21 and the fact r(—y) = r(y), there exist 
constants Cq > and < 5 < ah such that, 



dy 



1 


H-y) D ~ qI/a 


°dy< [ C 5 


r(y) D ~ qI/a 






Jr{y)<n 





<C e [ r(yr h ~^ S dy<^. 

Jt(«)<7 



r(j/)<7 

At the same time, (|2.1ip implies 

{y£R d : t(x + y) < 7} C {y G R d : r(y) < C ( 7 + r(-x))} 

= {yGM d :r(y)<Co(7 + T(x))}. 

Consequently we derive that 

/ ^x-y) -* 1 /* °dy= [ </>(- 

<C«[ T(y) ah - (1 ~ s dy<oo. 



' r{y)<C (l+r{x)) 

Combining the above shows that for any 7 > 

{x _ y) D-qI/a_ (j)( _ y) D-qI/ a 



r(j/)<7 



dy < 00. 



(4.9) 



Next we consider the integral on the set {y € M. d : r(y) > 7} for sufficiently large 
7 such that 4>{-y)- 1 T(x) < 1, Ci(/>(-?/)- /3 r(x) /3 < 1/2 and 0(-y) > 1. This is 
possible because of (|4.6p . Note that for any 3/2 > u > 1/2, from the fact 

dgg "" J/ " = ±eIn«(fl- 9 //a) = < D _ qI/a)s aD-(l +q /a)I 

ds ds 
and Lemma 14.21 there exists C7 > such that 



D — — 

a 



lVu 



lAu 



|| u Z>- ff //a _ j|| < 

<C 7 D 

Since 4> is E'-homogenous and 4>(—y) > 0, we have 

^_ y f- q I/a_^_ y) D-qI/ a 

H-y) D ~ qI/a 



, D -(i +9 / a )/,| ds 



(4.10) 



ti- II 



< 



4>{<jr\-y)x-<i>- E {-v)v) D - qI,a -i 



(4.11) 
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On the other hand, t(4> E (—y)x) = (f> 1 (—y)r(x) < 1 and <j>{—<j> E {—y)y) = 1, we 
can use (|4.7p and (|4.8p to derive 



; (- y )x - r s (-y)y) - 1 < Cj [r{4>- E {-y)x)\ 



(4.12) 



Since the last term is less than 1/2, we can apply (|4.10p with u = 4){4>~ E {— y)x — 
(j)~ E (-y)y). Hence, we derive from (|4.1ip . (|4.10[> . (|4.12|> and Lemma 1131 that for 
some < Si < (/3 - H)a 



<j>{p _ y) D- q I/a_^_ y) D- q I/ a 

<C? 4){-y) D - qI/a 

<c s H-y) D ~ qI/a 



a 


D-*I 






a 





\-y)x-<\r E {-y)y)-A 



■ 4>{-y)-^T{x)^ 



< C 9 cP(-y) aH+5l ~ q ~ aP T{x) al3 

< C 10 T{y) aH+8l - q - a(3 r(x) a(3 . 

This and Lemma |4. 1 1 yield 

f <p(x - y )D-*i/°> - 4>{-y) D - qI/a 

Vr(j/)>7 

< C u r{x) aP / r -<*(P-H)+h-i dr < ^ 



dy 



(4.13) 



Combining fj4.9f) and (|4.13p . we get (14. 5p which shows that X^, is well defined. 

(ii). To show the stochastic continuity of the a-stable random field X^, it is 
sufficient to verify that E(exp{z(0, X^(x + Xq) — X ( p(xo))}) — > 1 for all xq S M. d and 
9 G R m . By Theorem 2.2 it is enough to prove that for every xq £ R rf , we have 



dy as nO. (4.14) 



/ <Mz + x - ^-«'/« - 4>(x - y) D ~ qI / a 

By a change of variables, f|4. 14[) holds if 

T£(x) ^ as x -> 0. 

From the continuity of <^> and the continuity of the function (r, -D) -4 r D (see Meer- 
schaert and Scheflier [2DJ p. 30]), we have that 



(x - y) 



D-ql/a 



H-y) 



D-ql/a 



as X 



for every y £ M rf \{0}. Moreover, from the argument in Part (i), it follows that for 
a sufficiently large 7 > 0, there exists a positive constant C12 such that ||<£(x — 
y ^D- q I/a _ ^_ y ^D-qI/ a ^ ig bounded by 

*(*, y) = C 12 l {T(s/) < 7} (r(y) Q ^ + r(x - y) *"*"*) 
+C 12 T(y)^ +<5l -^-^r(x)^l {r( , )>7} , 
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where < 5 < ah and < Si < a(f3 — H). It is easy to see that $(x, y) — > $(y) a.e. 
as x — > 0, where 

$(y) = 2C 12 l {T{ y)< l} r(y) ah - 5 - q , 

and that 



/ $(z,y)dy 

JR d 



$(y)dy. 



By the generalized dominated convergence theorem (see [101 P-492]), (|4.14p holds, 
(iii). In order to show that for all r > 

{X^x), x€R d } = {r D X^(x), x G R d }, 

we note that, by Theorem 12.31 it is sufficient to prove that for all k > 1, Xj G R' 
and 0,- G R m (j = 1,2, • • • , fc) 



D-ql/a _ ^/_y\D-qI/a 
D-ql/a _ ^{ySD-ql/a 



dy 



(4.15) 



dy. 



This can be verified by an appropriate change of variables. By the .E-homogeneity 
of <b, we have 



sD-ql/a 



D-ql/a 



dy 



/ Y^L(r E X] -y) 

= / E f^~ 9//Q f^i - y) D_9//a - <K-y) D - 9//Q 

= / E - */) D -" /a - ^-y) D ~ qI/a 



dy 



r q dy 



dy. 



This proves (I4.15P and thus X is an (E, D)-o.s.s. random field. 

(iv). In the same way, we can verify that X${x) has stationary increments. The 
details are omitted. □ 



Finally, we prove Theorem 12.61 
Proof of Theorem \2.6l The proof is essentially an extension of the proofs of The- 
orem 4.1 and Corollary 4.2 in [3]. We only show that the stable random field X^ 
is well defined. Then properness of the X^ follows from the fact that for the ma- 
trix %l)(y)~ D ~ qI l a is invertible for every y G M d \{0}. The verification of the rest 
conclusions on X^ is left to the reader. 
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By Theorem 12.41 it suffices to show that 

T^):= / (\l-co S (x,y)r + \ S m{x,y)\ a ) m~ D ~ qI/a 



dy < oo. 



Let (ti(x), h(x)) be the polar coordinates of x under the operator E* . By (|2.12|) 
and Lemma 14.21 there exist < 5 < [j^(ai — H)] A (ah) and C13 > 0, such that 



and 



L {n (y)>l} 



l {n(y)<i} 



-D-ql/a 



< C 13 n(y) 



< Ci 3 ri(y) 



-ah+S—q 



-aH—8—q 



Hy)- D - qI/a 

Then by Lemma 14.14 Y i/j( x ) is bounded by 

° 13 !i I 1 ~ cos ^' r " B *^l° + |sm(x,r s ^)| a ) r ~ ah+& - V(d0)dr 
+C Vi 1^ J (jl - cos(x, r E *6) \ a + I sin(x, r E *9)\ a ^} r~ aH ~ s ~ 1 a(d9)dr. 

Note that there is a constant C14 > such that 

|l - cos(x, r E *6)\ a + I sm(x,r E '9)\ a < C 14 (l + \x\ a ) (r a{ai ~ S) A l). 
Therefore 



T£>(x) <C 15 (l + \x\ a )a(X 







/OO /*1 
r -a^+5-l dr+ / r a(ai-H)-(l+a)S-l dr 



Since < 5 < ~~ H)] A (a/i) and cr is a finite measure on So, we have 

T^(cc) < 00 for every x € ffi d . This proves that is a well-defined stable random 



field. 



□ 



The moving-average-type and harmonizable-type o.s.s. stable random fields are 
quite different (e.g., even in the special case of D = I, the regularity properties 
of and X^ are different.) From both theoretical and applied points of view, 
it is important to investigate the sample path regularity and fractal properties of 
the (E, D)-o.s.s. SaS'-random fields X^ and X^. We believe that many sample 
path properties such as Holder continuity and fractal dimensions of X^ and X^ are 
determined mostly by the real parts of the eigenvalues of E and D. It would be 
interesting to find out the precise connections. We refer to Mason and Xiao |19] , 
Bierme and Lacaux [3] and Xiao [32j [33] for related results in some special cases. 
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